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Abstract 

We construct an exact solution for the spherical gravitational collapse in a single coordinate 
patch. To describe the dynamics of collapse, we use a generalized form of the Painleve-Gullstrand 
coordinates in the Schwarzschild spacetime. The time coordinate of the form is the proper time of 
a free-falling observer so that we can describe the collapsing star not only outside but also inside 
the event horizon in a single coordinate patch. We show the both solutions corresponding to the 
gravitational collapse from infinity and from a finite radius. 

PACS numbers: 04.20.Jb, 04.40.Nr, 04.70.Bw, 97.60.L 
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I. INTRODUCTION 



One of the most important predictions of general relativity is black hole. Numerous 
studies have been done on the properties of the black hole and the formation by gravita- 
tional collapse. The theory of gravitational collapse was initiated in 1939 by Oppenheimer 
and Snyder as a spherical contraction model of a uniformly distributed dust star. The 
standard method in [1, is making a physically reasonable junction of the two different 
spacetimes corresponding to the interior and exterior regions of the collapsing body. The in- 
terior dust solution is given by the Friedmann- Robert son- Walker metric in the synchronous 
comoving coordinates. We impose a junction condition at the surface of the star so that the 
solution connects smoothly to the exterior Schwarzschild solution. The interior and exterior 
solutions are described in different coordinate systems. Although it is nothing wrong to con- 
struct solutions in such a manner, one cannot describe the dynamics of the collapsing star in 
terms of the coordinates of the observer outside the event horizon. The main purpose of this 
paper is to describe the both regions inside and outside the horizon by a single coordinate 
system in a physical way. 

The Painleve-Gullstrand coordinates used in this paper is, in fact, the key to a simple 
physical picture of black hole and gravitational collapse. It was first introduced by Painleve 

and Gullstrand [3] in 1921. We leave the details of its history to (sj, but mention here 
only its notable properties. Unlike the Schwarzschild form, the Painleve-Gullstrand metric 
tensor has an off-diagonal element so that it is regular at the Schwarzschild radius and has 
a singularity only at the origin of the spherical coordinates. In other words, the surfaces 
t = constant traverse the event horizon to reach the singularity. Therefore, the Painleve- 
Gullstrand coordinates are convenient for exploring the geometry of collapsing star and black 
hole both inside and outside the horizon altogether by a single coordinate patch. Moreover, 
the space given by the Painleve-Gullstrand coordinates can be intuitively regarded as a river 
whose speed of current is the Newtonian escape velocity at each point jsl. The velocity of 
the current hits the speed of light at the event horizon and then inside it there is nothing 
one can do to get sucked into the center of the hole. This feature results from the fact 
that this coordinate system adopts a time coordinate as measured by an observer who is at 
rest at infinity and freely falls straightforward to the origin. This physical picture provides 
a starting point of a method to describe gravitational collapse. In the present work, we 



generalize the Painleve-Gullstrand metric to incorporate gravitational collapse. 

The organization of the paper is as follows: in Sec. [TTl we introduce generalized Painleve- 
Gullstrand coordinates with the time coordinate being the proper time of a free-falling 
observer. In Sec. IIIII we solve the Einstein equation in the interior region of a spherical 
collapsing star of uniform perfect fluid, while in the exterior region the metric is given by 
the standard Painleve-Gullstrand form. In Sec. II VI we see the correspondence of the collapse 
from a finite radius and the model by Oppenheimer and Snyder. Sec. |V] is devoted to 
summary and discussion. 



II. PAINLEVE-GULLSTRAND COORDINATES 



In the case of spherical symmetry without charge, black holes are described by the 
Schwarzschild metric. The standard Schwarzschild coordinates (tg, r, 6, 0) are valid only 
outside the event horizon. However, the Painleve-Gullstrand coordinates {tp,r,6,(f)) enable 
us to describe also the inside of the horizon as well as the outside in a single patch of the 
coordinates. 

In this section, we derive the generalized Painleve-Gullstrand coordinates. We say here 
"generalized" in the sense that we introduce free-falling observers who start not only from 
infinity but also from other general points. The following mathematical derivation of the 
generalized Painleve-Gullstrand coordinates is close to the derivation in [gI] , but the physical 
situation and the time coordinate are different. 

The time coordinate tp of this family is the proper time r of an observer who freely falls 
radially from rest. Let us start with the Schwarzschild metric given in the standard form 

ds^ = -f{r)dU^ + f-\r)dr^ + {dO^ + sin^ Odcf?) , (1) 

where /(r) = 1 — 2M/r. In the Schwarzschild spacetime, the four- velocity u/" of an observer 
at the spacetime coordinates Xs^{t) is defined by Us^ = dxs^/dr = Xg^ with r being the 
proper time, which satisfies the normalization condition 

- 1 = g^uUsV = -fts + /-V'. (2) 

Note that the energy per unit rest mass of a test particle, 

e = -g,ui^u/ = fi, (3) 



is a constant of motion since = (d/dtg)^ is the timelike Killing field. Then the four- velocity 
Us^ can be explicitly written as 

ws'^ = (4, r, e, = (y , -V^'-f, 0, 0^ (4) 

and 



Us^ = gf.uUs" = { -e, J , 0, I . (5) 



Since we choose the Painleve-Gullstrand time coordinate tp as the proper time of the free- 
falling observer, the geodesic is orthogonal to the surfaces tp = constant and the geodesic 
tangent vector is equal to the gradient of tp: 

d_ 

that is to say, 



= tpi^sn, (6) 



_ f 

dtp = edtg + — — -j dr. (7) 

Consequently, the Painleve-Gullstrand metric takes the form 

ds^ = -dtp^ + ^{dr + v{r)dtpf + r^dVl^ (8) 

where 

vir) = - f{r) (9) 

is radially free-falling velocity. The observer in geodesic motion have the normalized four- 
velocity 

< = ftp, r, e, 0) = (1, -V, 0, 0) . (10) 



Note that the metric form ([8]) is different from that given by Martel and Poisson [61 , for our 
time coordinate is e times larger than theirs. This is because our metric is characterized by 
the free fall from various points at rest, while theirs by the free fall from infinity at various 
initial velocities. 

The important point to note is that the metric (|H]) indicates an analogue of the conser- 
vation of energy in the Newtonian mechanics, 

E=U^Y+^ir)^ (11) 



2 \dtpy 

where E = {e^ — l)/2 is a conserved energy and $(r) = —M/r is a gravitational potential 
energy of a central force field. In particular, if the particle freely falls from rest at infinity. 



the conserved energy E is zero (i.e., e = 1) and the metric ([8]) reduces to the standard form 
given by Painleve and Gullstrand: 
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ds^ = -dt/ + ^dr + y —dtpj + rdVL\ (12) 

For a free fall from infinity, the radial velocity v is the Newtonian escape velocity ^j2M/r. It 
is obvious from the metric ([H]) and ( fT2|) that the Painleve- Gullstrand coordinates are regular 
at the horizon r = 2M. This enables us to deal with the geometry of black hole both inside 
and outside the horizon. 

In the subsequent sections, we will consider the solution of the Einstein equation with 
matter. In general, the energy E and the mass M are functions of tp and r, not constant 
values. With the physical picture in mind and motivated by ([8]), we make the ansatz for the 
metric in the generalized Painleve-Gullstrand form 

ds^ = -dtp^ + ^^^^^^ (dr + v{tp, r)dtp) ' + r^dQ^ (13) 

where 



v{t„r) = ^2E{t„r) + ^-^^^^. (14) 

We may note, in passing, that our ansatz is not a particular case of^he metric given as 
a generalized form of the Painleve-Gullstrand metric by Lin and Soo ^] since their metric 
belongs to the category of the metric in 



III. SPHERICAL GRAVITATIONAL COLLAPSE FROM INFINITY 

In this section, we solve the Einstein equation in the spherical gravitational collapse. 
The standard form of the Painleve-Gullstrand metric gives a straightforward approach to 
the solution. 



A. Interior solution of the Einstein equation 

According to Birkhoff's theorem, the Schwarzschild solution is the only solution of the 
vacuum Einstein equation for a spherically symmetric spacetime. In particular, even if 
matter distribution is not static but moving in a spherically symmetric way, the exterior 

5 



vacuum region is given by the Schwarzschild metric. As shown in Sec. HTl we can express 
the Schwarzschild metric in the Painleve-GuUstrand form f lT2|) . Meanwhile, for the matter 
solution in the case of the gravitational collapse, the metric is assumed to be of the form 



ds' = -de +\^r+^j j + r''dn\ (15) 

that is, the generalized form f lT3|) with E = {]. In the following discussions, the Painleve- 
Gullstrand time coordinate is represented by t, instead of tp in Sec. [Tll Inside the spherical 
star, the Einstein equation, SnT^u = R^y — {1/2)6'^ ^R, reduces to 

8vrT% = (16a) 

8nT\ = (16b) 

^1 2m' 2m ( 2m\~^^'^ 

%'kT\ = + — — , 16c 

ji /™ l'2m\ mm' ( 2m^ 



8vrT^2 = 8vrn = -— +(£^ + -) ( — ) -^( — ) , (16d) 



m I m m 

7^ y V 7^ y ^2 \ qy 

where the dot represents the differentiation with respect to t and the prime the differentiation 
with respect to r. From equations (I16al) . (116bp and (I16cl) . there exists an identity, 

T^ = T% + nf— V^', (17) 



r / 

and then only three components of the Einstein equation are independent. 
For simplicity, we shall take the stress-energy tensor T^y of the perfect fluid: 

T^- = {p + P) u^u"" + Pg>"', (18) 

where p is the energy density, P is the pressure, and u'^ = (1, —v(t,r), 0, 0) is the four- 
velocity of the matter moving radially. Substituting equation ( !T8|) into equation (fT7|) . we 
obtain the radial velocity v{t,r) as the escape velocity 



2m(t,r) , . 

T;(t,r) = y^^. (19) 

Therefore, the four- velocity describes the radial free fall satisfying = (—1, 0, 0, 0) 
and u^Un = —1. Since it is clear from equation (fT9l) that the perfect fluid falls at the 
escape velocity, the physical picture is also clear — matter is attracted by the other matter 
distributed inside it. 
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Integration of equation fll6ap yields the mass function 

m{t, r) = An / p{t, r)r'^dr. (20) 
Jo 

Assume that the density has a separable form p = f{r)h{t) and that the pressure is given 
by P = (7 — l)p, where a constant 7 is usually greater than unity. Substitute them into 
equation (]16bp and impose the conditions rn\,^^Q = and pI^^q = 00, then the functions 
/ = 3A^/27r7^ and h = 1/9XH'^ are obtained, in which t < and A is an integration constant. 
Thus the matter density of the star is uniform: 

where t takes the value from —00 to 0. Furthermore, the other components of the Einstein 
equation reduce to a single equation 

P = SLZJ_, (22) 

It is consistent with the form that we supposed above. 

In the following subsections, we consider the simple case of dust (P = 0, i.e., 7 = 1) to 
explore the gravitational collapse in greater detail. 

B. Matching at the surface of the star 

In this subsection, the boundary surface between the interior dust region (r < R(t)) and 
the exterior vacuum region (r > R{t)) will be discussed. Let R{t) be the surface radius of 
the star at time t and M be the mass inside the surface. It is natural to impose a boundary 
condition at the surface r = R(t), that is, the mass, 

4:71 

m{t,r)l^^^^^ = —R^p = M, (23) 
of the star is constant, and then the radius of the boundary is given by 

Rit) = (^^Hry\ (24) 

This means that the surface of the star is at rest at infinity and its radius monotonically 
decreases to zero as t — )■ 0. In addition, the motion of the surface is geodesic. 



Now that the equation of the surface is determined, we are in a position to verify that 
the interior and exterior regions are smoothly matched each other. The interior solution 
(r < R{t)) given in Subsec. IIII Al has the line element 

ds't = -dt^ + [dr + v^{t,r)dtY + r^dfi^ (25) 

where the interior velocity 



is interpreted as the Hubble contraction flow ( — R{t) / R(t))r which is zero at the initial 
state and diverges at the final stage. At the surface r = R{t), this velocity becomes 

AM ^ 



and the interior line element becomes 

ds^^^^^^^ = -df + R\t)dQ^. (28) 
On the other hand, the geometry of the exterior region (r > R{t)) is 

dsl = -dt^ + {dr + v+{r)dtf + r^dn'^, (29) 
where the exterior velocity 

v,{r) ^ ^ (30) 

is the ordinary escape velocity in the Newtonian gravity. Hence when the radius r equals 
R(t), the interior and exterior velocities coincide with each other: v^{R{t)) = v^(t,R{t)) = 
V{t) and the line element of the exterior is the same as that of the interior: c?s^|^_^^^^ = 
ds'^\j._^^^y Moreover, the radial infall velocity V{t) of the surface of the star and the radius 



R{t) of it have a relationship 



4M ^ 



V{t) = ( ^ ) = (31) 



This implies that the surface motion is geodesic and the surface velocity hits the speed of 
light at time t = -AM/d,. 

Turning now to the matching of the extrinsic curvature of the boundary surface calculated 
both inside and outside the surface radius. Since the unit tangent vector to the geodesic 
curve of dust at the surface is 

«^ = (1, -y(t), 0, 0), (32) 
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the unit normal vector to the boundary surface is given by 

= {V{t), 1, 0, 0) (33) 

which satisfies the orthogonahty w'^ra^ = 0. Using the normal vector ra^, the extrinsic 
curvature is 

K,, = -h,"K^Vpn^, (34) 

where h^i, = g^i, — n^rii, is the projection to the surface. Consequently, the components of 
the extrinsic curvature of the boundary surface are 

K22 = = -R{t), (35a) 

K^, = Kt3 = -R{t) sin' 9, (35b) 
and all the others components vanish, especially 

K^o = -V+ v±l=ni,) - ^ ^±L=iiW = 0- (35c) 

Therefore, the interior region and the exterior region are smoothly matched at the boundary 
surface r = R{t), and then we can describe the geometry of all the spacetime by a single 
coordinate system {t,r,9,(f)). 

C. Causal structure 

The previous Subsec. UlI BI focuses on the local property, i.e., the smoothness of the metrics 
at the boundary surface. In this subsection, we study the global structure of spacetime with 
the collapsing dust considered above. The Painleve-Gullstrand coordinates manifestly give 
a physical picture of the whole spacetime of the gravitational collapse. 

The interior solution (125|) expresses a contracting fiat Friedmann universe. The dust 
region is inside the boundary surface r = R{t) that freely and radially falls at the speed 
V{t) with time —00 < t < 0. In this region, there exists the event horizon and the apparent 
horizon defined by r = 3t + 3R{t) and r = —3t/2, respectively. The event horizon emerges 
from the origin r = at the time t = — 9M/2 and its radius becomes larger until it equals 
to the surface radius R{t). The time that the two radii become equal to each other is 
t = — 4M/3 and the size is r = 2M, as expected. After the surface of the star is hidden by 
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FIG. 1. The Penrose diagram for the collapsing dust star. The shaded area is the collapsing 
dust. The thick curve represents the trajectory of the boundary surface r = R{t). The dotted and 
the dashed-dotted lines represent the apparent horizon, and the event horizon, respectively. The 
trapped region is surrounded by the apparent horizon. 

the event horizon, the timelike apparent horizon appears at the same intersection. Outside 
the apparent horizon r > — 3t/2, surfaces r = constant are not timelike but spacelike, and 
the infall velocity exceeds the speed of light V-{t,r) > 1 so that this region is a trapped 
region. And the contracting dust star collapses into a singularity at t = 0, consistent with 
the singularity theorem. 

The whole spacetime diagram for the gravitational collapse is given by the interior dust 
region mentioned above and the exterior Schwarzschild region in FIG. [U It expresses the 
dust collapsing from rest at infinity and the eventual black hole. The dust fiat Friedmann 
universe case (P = 0) for the interior spacetime was investigated in in the Painleve- 
Gullstrand coordinates. Here we have studied the same problem from the new viewpoint of 
the generalized Painleve-GuUstrand metric in which physical picture is clear. 



IV. SPHERICAL GRAVITATIONAL COLLAPSE FROM A FINITE RADIUS 

In the previous section, we have discussed the description of the spherical gravitational 
collapse that starts from infinity. Its interior metric (|25|1 is expressed by the fiat Friedmann 
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universe. However, the standard model of a collapsing star, including the one considered by 
Oppenheimer and Snyder is collapse from rest at a finite initial radius. It is not trivial 
to apply the idea of the Painleve-GuUstrand coordinates to the situation of the collapse 
starting with a finite radius. The standard metric of this type is that of a contracting closed 
Friedmann universe as a solution of an interior dust. Although the similar solution is 
given by Gautreau and Cohen it is complicated even in the vacuum region since it is 
given by implicit functions of time and radial coordinates. 

In this section, we briefiy show the solution of the collapse from a finite radius as an- 
other particular case of the generalized Painleve-Gullstrand form flTS]) . To begin with, we 
consider the boundary surface that freely falls from rest at a radius Rq. At the surface, the 
conservation law f lTT]) gives the energy 



M , , 

^ = -^ (36) 



and the infall velocity 



where Rify is the surface radius smaller than the initial radius -Rq- 

Next, when we assume that the exterior solution ranges from the contracting surface 
radius i?(t) to the initial radius i?o (i-e., -R(t) < r < .Ro), the energy remains the same as 

iSD, 

^+ = (38) 
and therefore the infall velocity becomes 



/2M ^ /2M 2M , , 

in the exterior region. The exterior metric is therefore given by 

ds\ = -dt^^^ {dr^Ji—-'^-^dt\ +r'dQ'. (40) 

2mv V ^ ^0 y 



Finally, in the interior region < r < R{t), since the energy is now a function of time 
and radius coordinates, we make the ansatz 

m(t, r) 

E4t,r) = (41) 

r 
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for the energy and 



„_(Vr)^;^.2^_,,,-)^^^(l-^) (42) 

for the velocity in the generahzed Painleve-Gullstrand metric f ll3p . Note that (155]) and 
E^{t,r) f HT]) coincide with E f l5B]) . and also that f+(r) ( 15^ and f_(t,r) coincide with 
(137]) at the surface r = 
In the case of the uniformly distributed dust, we can solve the Einstein equation taking 
into account the boundary condition which is the same as that in Sec. IIII B[ The mass 
function reduces to 

m{t,r) = —r'p{t), (43) 
3M 



where the surface radius is 



the time coordinate is 



R(t) = ^{l + cosv), (45) 



^ = V^(^ + sin^) (46) 



and the parameter rj takes the value from to vr. This interior solution 



dsi = -dt' + ^ ( dr + - —-fTTT^dt 1 + r'dn' (47) 

2M ( r ^ 

1 - 



2M 


2M 


R{t) 


Ro 



2 



Ro \Rit) 

expresses the contracting closed Friedmann universe. As before, one can verify that the two 
metrics corresponding to the interior and exterior solutions smoothly match at the boundary, 
by using the normal vector 

^/. = ^7^i^(n^), 1,0,0). (48) 
Namely, the components of the extrinsic curvature of the boundary surface are 



^22 = ^2^2 = -Vl + 2ER{t), (49a) 



7^33 = 7^3+3 = -Vl + 2ER{t) sin^ 9, (49b) 
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and all the others vanish. The interior velocity (H2|) is also characterized as the Hubble 
contraction flow ( — R(t)/R{t)^r, just as the velocity f l26|) in the Hubble flow in the flat 
Friedmann universe. In fact, when the initial radius is large enough to satisfy Rq ^ M, the 
energy and the velocity in the both regions become the same as those of the collapse from 
inflnity. 

In the further exterior region r > Rq, the metric is simply given by the standard 
Schwarzschild form up to the scale of the time coordinate: 

2M 

1 ^^2 

dsl^ = ^dt^ + + r^dn^ (50) 

2M 2M 
1 _ 1 

Ro r 

which smoothly matches with the exterior solution f HOj) at the boundary r = Rq. Note that 
the infall velocity f+(r) (15^ in the exterior region vanishes at r = Rq. Similar expression 
can be found in [lo| with a = - 2M/r)/(l - 2M/Rq) and E = -2M/r in the notation 
there. Rigorously speaking, we need another coordinate patch for r > Rq, giving up a 
single coordinate system in the entire region of spacetime. However, we are interested in the 
coordinate system which describes the inside and outside of the horizon, but not very much 
in the region r > Rq. 




t = o" 



FIG. 2. The Penrose diagram for the collapsing dust star of the initial radius Rq. The shaded area 
is the collapsing dust. The dashed line represents the surface t = as the initial hypersurface. 
The thick curve represents the trajectory of the boundary surface r = R{t). The dotted and the 
dashed-dotted lines represent an apparent horizon and an event horizon, respectively. The trapped 
region is surrounded by the apparent horizon. 
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The Penrose diagram for the dust star collapsing from the initial radius r = Rq is shown 
in FIG. m The surface radius r = R{t) is at rest at the initial time t = 0, i.e., rj = 0, 
and monotonically decreases to zero as 77 — t- tt. The event horizon emerges from the origin 
at ?7 = c, where c is a constant and equals to arccos(4M/i?o ~ 1) ~ arcsin ^2M/Ro, and 
subsequently its radius grows as r = a/ Ro/2MR{t) sm{ri — c). At rj = c + arcsin y^2M/RQ, 
the radius of the horizon becomes equal to the surface radius, 2M. After that, the event 
horizon (the outer apparent horizon) appears at r = 2M and the (inner) apparent horizon 
at r = ^R^{t)/2M. The regions of the infall velocity f+(r) > a/1 + 2i?+ and V-{t,r) > 
^/l + 2E_{t, r) are the trapped regions. 

So far is the metric for t > 0. Although we may continue the spacetime backward to the 
region t < sustaining the radius of the star Rq, this is out of the scope of the present work. 
Here, we simply to give the initial data for the collapse situation: at the hypersurface t = 0, 
take the uniformly distributed spherical dust of radius r = Rq and the exterior Schwarzschild 
region for r > Rq, as indicated in FIG. |2] 



V. SUMMARY AND DISCUSSION 

For the description of gravitational collapse of a dust star, we have introduced the gen- 
eralized Painleve-Gullstrand coordinates with the time coordinate being the proper time of 
a free-falling observer. We gave the solutions of the Einstein equation in the cases of the 
collapse from a finite radius as well as from infinity. The metric describes both the interior 
and exterior regions of the star, which smoothly match at the surface of the star. More pre- 
cisely, the metric is of class, while the metric component is of class. The choice of 
the Painleve-Gullstrand time coordinate enables us to write the solutions inside and outside 
the event horizon in a single coordinate patch. 

As another application of the Painleve-Gullstrand coordinates, we can consider a spherical 
model of expanding universe which consists of regions of void and inhomogeneous dust 
distribution. In a sense this is the inside-out version of the gravitational collapse in the 
present paper. 

14 



ACKNOWLEDGMENTS 



The authors are supported by Global Center of Excellence Program "Nanoscience and 
Quantum Physics" at Tokyo Institute of Technology. 



[1] J. R. Oppenheimer and H. Snyder, Phys. Rev. 56, 455 (1939). 

[2] C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravitation (Freeman, San Francisco, 1973). 

[3] P. Painleve, C. R. Acad. Sci. (Paris) 173, 677 (1921). 

[4] A. Gullstrand, Arkiv. Mat. Astron. Fys. 16, 1 (1922). 

[5] A. J. S. Hamilton and J. P. Lisle, Am. J. Phys. 76, 519 (2008). 

[6] K. Martel and E. Poisson, Am. J. Phys. 69, 476 (2001). 

[7] C.-Y. Lin and C. Soo, Phys. Lett. B 671, 493 (2009). 

[8] R. J. Alder, J.D. Bjorken, P. Chen and J. S. Liu, Am. J. Phys. 73, 1148 (2005). 
[9] R. Gautreau and J. M. Cohen, Am. J. Phys. 63, 991 (1995). 
[10] P. D. Lasky and A. W. C. Lun, Phys. Rev. D74, 084013 (2006). 



15 



